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I. INTRODUCTION 

In the last decade, the interest for ground state nu- 
clear properties, in particular, static moments, was re- 
newed due to modern Radioactive Ion Beam facilities 
which provide access to long chains of isotopes, including 
the radioactive ones, in their ground and isomeric states. 
Such nuclei distant from the /3-decay stability valley and 
often very close to the drip lines are of great interest for 
nuclear astrophysics. Spectroscopy techniques using high 
intensity lasers allow for precision measurements of nu- 
clear spins, magnetic and quadrupole moments which are 
much more delicate nuclear characteristics than masses 
and charge radii. As a result, the bulk of the data on nu- 
clear static moments becomes very extensive and compre- 
hensive [Ij creating a challenge to nuclear theory. Recent 
precise measurements of magnetic and quadrupole mo- 
ments of a long chain of copper isotopes from iV = 28 
to iV = 46 presented and analyzed in the comprehen- 
sive article [4I is a new evidence of this interest. For 
nuclear magnetic moments, this challenge was partially 
responded recently 0, |4|. Data on magnetic moments 
of more than one hundred odd nuclei were reproduced 
with an accuracy of 0.1— 0.2/iN within the self-consistent 
Theory of Finite Fermi Systems (TFFS). In this study, 
the QRPA-like TFFS equations for the nuclear response 
to the external magnetic field were solved on the base 
of the Energy Density Functional (EDF) by Fayans et 
al. 043- Especially high accuracy was reached for semi- 
magic nuclei considered in the "single-quasiparticle ap- 
proximation" where one quasiparticle in the fixed state 
A = (n, I, j, m) with the energy ex is added to the even- 



even core. Such nuclei where only one kind of nucle- 
ons is superfluid rule, spherical. In addition, 

the lowest 2+-state has usually an excitation energy of 
the order 1 MeV, the corresponding collectivity is not 
extremely high and the particle-phonon coupling is not 
sufficiently strong to destroy the single-quasiparticle ap- 
proximation. According to the TFFS Q, a quasiparticle 
differs from a particle of the single-particle model in two 
respects. First, it possesses a local charge and, second, 
the core is polarized due to the interaction between the 
particle and the core nucleons via the Landau-Migdal 
(LM) amplitude. In other words, the quasiparticle pos- 
sesses the effective charge Ceff caused by the polarizability 
of the core which is found by solving the TFFS equa- 
tions. In the many-particle Shell Model 9], a similar 
quantity is introduced as a phenomenological parameter 
which describes polarizability of the core consisting of 
outside nucleons. For nuclei with both non-magic sub- 
systems corrections to the single-quasiparticle approxi- 
mation due to 2+-phonon coupling are important and 
should be taken into account to describe magnetic mo- 
ments correctly [1, 0] • 

For quadrupole moments, the systematic calculations 
we know concern only the medium atomic weight nuclei 
of the 2p - If shell [| and 2p, I/5/2, Igg/z shell 0. Evi- 
dently, the first step in this direction for heavy nuclei was 
made in our recent article flo| where we analyzed effects 
of the density dependence by studying the 2^-excitations 
in two isotopic chains of semi-magic even-even tin and 
lead nuclei within a version of the self-consistent TFFS 
which is very similar to that in 0, H- The only difference 
is the use of a new version DF3-a |Tl| of the original DF3 
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functional [1, 0| • It differs from DF3 only by the spin- 
orbit parameters >f, V and the first harmonics gi , g[ of 
the spin LM amplitude. The corresponding force compo- 
nents determine mainly the spin-orbit doublet splitting 
and are of especial importance for high j-levels. The last 
two terms modify the spin-orbit component of the mean 
field together with the spin-orbit density Psiir) which 
is determined with particles of partially filled spin-orbit 
doublets. As it is well known, the first spin harmonics 
contribute to the mean field in a combination with tensor 
forces, see e.g. [lH. Therefore, they should be considered 
as the effective first harmonic or, on equal footing, the ef- 
fective tensor force. It should be mentioned that in both 
the functionals DF3 and DF3-a the parameter gi is taken 
equal to zero, thus, the spin-orbit splitting is determined 
with the set of three parameters, }i,>c',g'i. For brevity, 
we will name them just spin-orbit. The functional DF3-a 
is characterized by a rather strong effective tensor force. 

In this paper we concentrate on quadrupole moments. 
We limit ourselves to semi-magic nuclei, where the one- 
quasiparticle approximation is expected to provide a 
good accuracy. In addition to two isotopic chains con- 
sidered in [lO|, we calculate quadrupole moments of odd 
neighbors of even isotones with iV = 50 and iV = 82 and 
of four odd neighbors of the magic nucleus ^°Ca. We 
compare predictions of two functionals, DF3 and DF3-a, 
which have different spin-orbit components. It is impor- 
tant for the problem under consideration as quadrupole 
moments Qx are proportional to the Bogolyubov factor 
fl|: 



ul-vl = {ex - ti)/Ex, 



(1) 



/i)2-|-A|, with obvious notation. 



where Ex = ^JJel 
In the quadrupole moments problem we deal with the 
ground state of an odd nucleus or with very low-lying 
excited state when often the inequality \ex — P-\ ^ A^ is 
true. In such a situation, the quadrupole moment value is 
very sensitive to the single-particle energy ex- In its turn, 
the latter is very sensitive to the spin-orbit parameters 
of the EDF. This could help in fixing these parameters 
not known sufficiently well up to now. The reason is that 
nuclear masses and radii used mainly for finding effective 
force parameters are not-sensitive to such parameters as 
'^'j 51)51- Note that the relevance of the spin-orbit and 
effective tensor force to other low-ener gy n uclear phe- 
nomena was discussed recently in Refs. |14l4l6| . In this 
paper, we consider only the surface kind of pairing as 
motivated by previous research of Refs. 0, [TgI [l7j. 



II. BRIEF CALCULATION SCHEME 

The calculation scheme of the self-consistent TFFS 
based on the EDF method by Fayans et al. is described in 
detail in Ref. [lo| . Here we summarize several formulas 
which are necessary for understanding the main ingredi- 
ents of the approach. The EDF method by Fayans et 
al. [^-01 is a generalization for superfluid systems of the 



original Kohn-Sham EDF method [18|. For condensed 
matter case, similar generalization was carried out in Ref. 
[l9l |. In this method, the ground state energy of a nucleus 
is considered as a functional of normal and anomalous 
densities, 



£[pn (r) , pp (r) , Vn (r) , Vp {v)]d^r. 



(2) 



Within the TFFS, the static quadrupole moment Qx 
of an odd nucleus with the odd nucleon in the state A 
can be found in terms of the diagonal matrix element 
{\\V{ijj = 0)|A} of the effective field V in the static ex- 
ternal field Vb = y^l6n/5r^Y2o- In systems with pairing 
correlations, equation for the effective field can be writ- 
ten in a compact form as 



(3) 



where Cq is the local quasiparticle charge with respect to 
the external field Vq and all other terms are matrices. In 
the standard TFFS notation Q, we have: 



V = 






J" = 



C{io) Mi{uj) M2iuj) 



(4) 



(5) 



(6) 



where C, A4i, and so on stand for integrals over e of the 
products of different combinations of the Green function 
G(e) and two Gor'kov functios F^^) (e) and F^^) (e) . They 
can be found in Q. In our case, the local charges in Eq. 
© are = 1, = 0. 

Isotopic indices in Eqs. (jHH]) are omitted for brevity. 
The explicit form of the above equations is written down 
for the case of the electric (t-even) symmetry we deal 
with. In Eq. ([5]), J-" is the usual LM amplitude. 



J": 



Sp^' 



(7) 



is the density-dependent effective pairing interaction, 

(8) 



and the amplitudes = T^'^ stand for the mixed sec- 
ond derivatives. 



SpSiy 



(9) 
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In the case of volume pairing, one has F'^^ = 0, whereas 
for the case of surface pairing we deal the amplitude J^"^ 
is non-zero and should be taken into account when Eqs. 
dHH]) are solved. As the analysis of Ref. shows, 
the component V of the vector V', as a rule, dominates. 
However, the fields di,d2 also contribute, and sometimes 
significantly, to the value of Qx. 

The final expression for the quadrupole moment of an 
odd nucleus is as follows [l^l 

QT = {nl-vl)vr, (10) 

where v\ are the Bogolyubov coefficients and 

Vx = / V{r)Rl,^{r)rHr. (11) 

It should be noted that we use here, just as in Ref. [lOj . 
the diagonal in A approximation for the gap, A^.a' = 
Aa(5a,a' where the Bogolyubov set of equations is reduced 
to the BCS scheme. In the general case, the relation for 
the quadrupole moment similar to Eq. ()10|) is more com- 
plicated. The j-dependent factor in pip appears due 
to the angular integral [20^1. For j > 1/2 it is always 
negative. For odd neighbors of a magic nucleus the "Bo- 
golyubov" factor in PH)) reduces to 1 for a particle state 
and to —1 for a hole one. If the odd nucleon belongs to 
the superfluid component, the factor (ua~^a) ™ ^'i- fTO|) 
becomes non-trivial, see Eq. ^ . It changes permanently 
depending on the state A and the nucleus under consid- 
eration. Note that in the case of magnetic moments the 
factor of {u\ + v1) = 1 appears [ij in the relation analo- 
gous to (fTO|) . In our case, this factor determines the sign 
of the quadrupole moment. It depends essentially on val- 
ues of the single-particle basis energies £a reckoned from 
the chemical potential fj, as. Keeping in mind such sen- 
sitivity, we found this quantity for a given odd nucleus 
{Z,N + 1) or {Z + 1,N), N,Z even, with taking into 
account the blocking effect in the pairing problem [l3| 
putting the odd nucleon to the state A under considera- 
tion. For the V\ value in Eq. pil)) for superfluid nuclei, 
we used the half-sum of these values in two neighboring 
even nuclei. 



III. CALCULATION RESULTS 

The set of nuclei analyzed includes, in addition to odd 
neighbors of two chains of even is otop es of tin {Z ~ 50) 
and lead {Z = 82) considered in [MM: odd neighbors of 
even isotones with iV = 50 and = 82. To check the 
method for lighter nuclei, we calculated also quadrupole 
moments of odd neighbors of the magic ^°Ca nucleus. We 
concentrate on one of the two versions of pairing interac- 
tion considered in [l3|, the surface one, as it was favored 
by that investigation. On the other hand, we consider 
two versions of the EDF. In addition to the DF3-a func- 
tional used in fic| we made alternative calculations with 



the original DF3 functional of [d, m . All parameters of 
the two functionals are given in [lO|. 

Let us begin by analysing nuclei with odd neutron 
numbers. Quadrupole moments of the odd Sn isotopes 
are displayed in Fig|TJ Predictions for the two versions of 
the functional, DF3 and DF3-a with different spin-orbit 
parameters, are compared with each other and with the 
experimental data. In average, the results of both ver- 
sions reasonably agree with the data. Each theoretical 
curve crosses the zero line due to vanishing of the cor- 
responding Bogolyubov factor in Eq. (fTU)) . The anal- 
ogous experimental curve could be approximately dis- 
played only for the ll/2~-state. We see that it crosses 
the zero close to both theoretical curves, the difference 
between those is 5 A ~ 2. 

Before analysing the results for both versions of 
the EDF in detail, it is instructive to discuss the A- 
dependence of the position of the state l/iii/2 and com- 
pare it with the experimental one. For the tin chain, this 
"intruder" state plays a special role due to the high j 
value and closeness to the Fermi level. Indeed, the prop- 
agator £{uj = 0) of Eq. ^ playing the main role in the 
TFFS equation for the effective field V2+ {oj — 0) contains 
diagonal elements C\\ = — A|/(2i?|). Their contribu- 
tion is proportional to the {2j\ + 1) factor. If the energy 
E\{hiif2) is small, the contribution of this level to £ can 
be rather big. Therefore a change in the position of this 
level may change the result for the quadrupole moment 
significantly. 

As it was explained above, the DF3 and DF3-a func- 
tionals differ with spin-orbit parameters only, the dif- 
ference being important mainly for high j-levels. The 
1^11/2 level position accounted from the ground state in 
Sn isotopes, Ex — fi, for both the functionals for even 
tin isotopes is displayed in Fig. [5] in comparison with the 
experimental excitation spectrum in neighboring odd nu- 
clei. The corresponding values of the Bogolyubov factor 
u\ — Vx are shown in the right-up sector of the figure. 
Dealing with an even-even nucleus, say, even '^Sn 
isotope, there is a dilemma how to interpret a state |A) 
under consideration, either the hole state (i.e. the exci- 
tation of the '^^^Sn nucleus) or the particle one (the ex- 
citation of the ^"""^Sn nucleus). We use a simple ansatz: 
the state \i) is considered as a hole state if the inequality 
{u\ — vf) < takes place and as a particle state other- 
wise. Note that in the case of (m| — vf) — the differ- 
ence between the particle and hole energies is, as a rule, 
quite small. Fortunately, for the main part of the chain 
there is no contradiction between both the functionals in 
the "particle/hole" interpretation of the lft.11/2 state. In 
the left part, up to ^^'^Sn nucleus, we deal with particle 
states and, correspondingly, for each A value here the 
green triangle shows the excitation energy of the (A+l)- 
th isotope. On the right side, A > 120 experimental data 
show the ground state ll/2~, or almost ground state in 
I30,i32gj-^ isotopes. Therefore, a formal contradiction in 
this point for the ^^^Sn doesn't work. We see that both 
functionals predict ^-dependence of the excitation en- 



4 




FIG. 1: (Color online) Quadrupole moments of odd tin isotopes with DF3 and DF3-a functionals. Triangles with the bars 
indicate experimental data JJ. 
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FIG. 2: (Color online) For the /111/2 level in Sn isotopes: 
the energy E\ accounted from the ground state(the left-down 
legend; left axis) and values of the Bogolyubov {u\—v'l) factor 
(right-up legend; right axis). 



FIG. 3: (Color online) Single-particle spectrum of the ^^*Sn 
nucleus. Value of the Bogolyubov (it J — ) factor is given 
near each theoretical crossbar. 



ergy of the ft.11/2 level qualitatively similar to the exper- 
imental one, falling in the value with A increasing, but 
both theoretical curves are higher than the experimental 



one. Beginning from A — 124, the DF3-a points prac- 
tically coincide with the experimental ones. The DF3-a 
curve is above the DF3 curve to the left of the crossing 
point between A = 120 and A = 122, to the right, the 
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TABLE I: Characteristics of single-particle states in "»Sn 
nucleus. 



FIG. 4: (Color online) Effective field V{r) in "**Sn nucleus 
with DF3 and DF3-a functionals. 



position is inverse. For almost all the A values the in- 
equality {uj^—vf)BF3 < (w^ — w^)DF3-a is valid. Absolute 
values of the Bogolyubov factor for these functionals dif- 
fer often significantly, as a rule, more than the quadrupole 
moment values, see Fig. [T]and Table HIl Analysis shows 
that it occurs due to the opposite effect in the matrix 
elements V\ of the effective field. 

Let us analyze this point in more detail for the nucleus 
^^^Sn which is in the middle of the tin chain. The exci- 
tation spectrum for this nucleus is displayed in Fig. [31 
According to the ansatz formulated above, position of ex- 
perimental levels for 3/2+, 11/2^ states are taken from 
the spectrum of the ^^^Sn nucleus, that of 7/2+, 5/2+ 
states, from ^^^Sn isotope. Values of the Bogolyubov 
m| — v"^ are also calculated in corresponding odd nuclei 
with the blocking effect taken into account as described 
in Chapter 2. In this nucleus positions of 11/2^ state 
for both functionals are rather close to each other as this 
nucleus is in the vicinity of the crossing point in Fig. [21 
However energies of the 7/2+ state, also with rather big 
momentum value, are strongly different. On the whole, 
the density of states at the Fermi surface for the DF3 
functional is higher than for DF3-a one. This difference 
is, in fact, significant if one takes into account the factor 
of (2j -I- 1) with which each j-state comes to this quan- 
tity. To avoid misleading, we note that the comparison 
of the theoretical spectrum in the even-even A-nucleus 
which is the basis for the QRPA-like calculation with 
the experimental excitation spectra of neighboring odd 
nuclei, as is made in Figs. 12131 is, in general, rather ap- 
proximate operation. In heavy nuclei, corrections to the 
mean field single-particle spectra appear due to particle- 
phonon coupling. Phenomenological functionals, as DF3 
or DF3-a, or any SHF functional, take into account such 
corrections only on average. For quantitative description 
of single-particle spectra, all fiuctuations, from one nu- 
cleus to another and from one state to another, should 
be calculated explicitly. 
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FIG. 5: (Color online) Quadrupole moments of odd lead 
isotopes with DF3 and DF3-a functionals. Triangles with the 
bars indicate experimental data [l[. 



To complete the analysis of quadrupole moments of the 
odd neighbors of ^^®Sn, we presented in Table [J matrix 
elements V\ entering Eq. (jlOp and quadrupole moments 



Q\ = (u| — vj^)V\. Note that there is a small difference 
between these values of Q\ and those in Fig. [Hand Table 
im as in our systematic calculations of quadrupole mo- 
ments we used in this product for an odd nucleus under 
consideration the half-sum of V\ values in neighboring 
even nuclei. We see that absolute values of V\ for the DF3 
functional are almost two times bigger than for the DF3-a 
functional used by us in the previous article fld\. How- 
ever, in three cases of four the difference of Bogolyubov 
factors is opposite and compensates significantly too big 
values of V\ for the DF3 functional. This compensa- 
tion is almost complete for 7/2+ and 11 /2~ states, the 
corresponding quadrupole moments for two functionals 
almost coincide in vicinity of ^^*Sn. This is not the case 
for the 5/2+ state where two Bogolyubov factors almost 
equal to each other and yx[DF3] almost twice bigger of 
yx[DF3-a)]. A strong enhancement of the effective field 
V{r) for the DF3 functional in comparison with DF3-a 
one is shown in Fig. [4l It is seen that the surface maxima 
of Vp{r) and Vn{r) functions are almost twice higher in 
the DF3 case. The higher density of states at the Fermi 
surface discussed above is causing this effect. 

For a deeper understanding of this enhancement effect. 
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FIG. 6: (Color online) Quadrupole moments of odd-neutron 
neighbors of even A*" = 50 isotones with DF3 and DF3-a func- 
tionals. Triangles with bars indicate experimental data jj. 



we calculated for the nucleus ^^^Sn characteristics of the 
2^ state with the DF3 functional, uj2 = 0.913 MeV and 
B{E2) — 0.272 e b, jmd compare them with those cal- 
culated previously (lOl for the DF3-a functional, uj2 = 
1.217 MeV and B{E2) = 0.172 e b. Note that the DF3- 
a predictions are in reasonable agreement with experi- 
mental data, thus the DF3 functional leads to too col- 
lective 2^ state. It explains qualitatively an additional 
enhancement of the static quadrupole field for the DF3 
functional. Indeed, at small external quadrupole field fre- 
quency w, the pole in the uj — UJ2 point term dominates 
in the effective field. 



^(^) = — 



(12) 



where the regular term Vr{lu) ~ const at small w. The 
residue of the pole term is Xs = ffs being the 

amplitude of creating the s-th {2^ in our case) phonon 
Q. It can be easily shown that the relation X2 oc B{E2) 
is valid. If to neglect for an estimation the regular term 
in Eq. (HI]), we arrive to V{uj = 0) oc B{E2)/u}2- After 
substituting to this relation the values given above, we 
find this ratio equal to 0.298 for the DF3 functional and 
0.141 for DF3-a. This explains the effect under discus- 
sion. 

To evaluate the agreement with experiment quantita- 
tively, we calculate the mean theory-experiment differ- 
ence 



{SQ) 



2 

rms 



(13) 



with obvious notation. F or the tin chain under considera- 
tion, Af = 11, we have a/ (5Q)r, 



0.167 eb for the DF3 



functional and y ((5Q)rins = 0.164 e b for the DF3-a one. 
Thus, for the DF3-a functional agreement is somewhat 



better. A remark should be made concerning practical 
application of Eq. p^ . All terms of the sum come with 
the same weight independent on the experimental error 
of each measurement. The reason is that, as calculations 
have shown, our theoretical accuracy for quadrupole 
moments is not better than 0.1 e b. Practically all the 
data we use for comparison have better accuracy there- 
fore for our analysis they are equivalent. The nucleus 
209pj^ is the only exception where the experimental error 
is 0.2 e b. To take it into account, we multiply this term 
with the factor (0.1/0.2)^ = 1/4. If the table |l| contains 
several different data on the quadrupole moment under 
consideration, say 3, we take into account all of them in 
the sum of (fT5|) with the weight 1/3. 

Let us go now to odd lead nuclei. Fig. [5j We see that 
again there is a reasonable agreement with experiment for 
both the functionals with the exception of the intruder 
state lii3/2- For the latter, in the case of the DF3-a func- 
tional, a reasonable agreement with the data is found 
only for two lightest isotopes i9i'i93p^ p^j. Q^^j^gj- fom- 
nuclei with known values oiQ{lii^^2) disagreement is sig- 
nificant, the maximal value of 6Q = Q"^ — Q°'^p for ^^^Pb 
is 0.60 e b. For the DF3 functional the situation is notice- 
ably worse. Now the rms deviation of the theory from ex- 
periment, Eq. (Hg, is \/((5Q)2„3[DF3-a] = 0.277 e b and 



(^Q)2^,[DF3] = 0.386 e b. Now the quality of agree- 
ment became significantly worse than for tin chain for 
both functionals but the deterioration is much stronger 
for the DF3 functional. 

In Fig. |6] quadrupole moments are displayed of odd- 
neutron neighbors, N = 50 ± 1, of even isotones with 
N = 50. In this case, the proton-subsystem is superfluid 
and the neutron Bogolyubov factor in Eq. ((TU)) is ±1. 
Agreement with the data is almost perfect for the DF3- 

a functional, y ((5Q)j2,„s[DF3 — a] = 0.041 e b, and only 



a little worse for the DF3 functional, y ((5Q)rj,js[DF3] = 
0.116 e b. For this chain we showed also quadrupole 
moments of the intruder state I/111/2 although none of 
them was measured. Contrary to tin and lead isotopic 
chains, none of the functions Q\{A) changes the sign. 
The situation is similar for odd-neutron neighbors of even 
isotones with iV = 82 presented in Fig. [7l Again agree- 
ment with the data is very good for the DF3-a functional, 

DF3 



a] — 0.093 e b, and again only a little 



worse for the DF3 functional, ((5Q)2^,JDF3] = 0.119 e 
b. In this case, there is one experimental data for the 
intruder state lii3/2 for which the DF3- value practically 
coincides with experiment. 

Table |ll] contains odd-neutron nuclei with known ex- 
perimental quadrupole moments. We have seen that 
for each chain considered the DF3-a functional describes 
data better than the DF3 one. Therefore we present here 
only predictions of the DF3-a functional. The table is 
divided in two parts. In the bottom, the intruder states 
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1^11/2 and lii3/2 are collected. The rms deviation of 
the theoretical predictions from the data calculated with 
Eq. [13] for all 42 moments in this Table is rather big, 

\/('5Q)?ms[DF3 - a] = 0.189 e b. For the DF3 functional, 

the error is bigger, •y/(5Q)^[DF3] = 0.240 e b. Examin- 
ing the table, we find that the main part of big deviations 
SQ is concentrated in the "intruder" part of the Table 
containing 15 moments. If to make the calculation 

only for the latter, we find {SQ)^^^g [intruder] — 0.269 e 
b. At last, rather moderate value of the error will occur 
if we limit ourselves with quadrupole moments of the rest 

of 27 "normal" states, (SQ)^^^ [normal] = 0.125 e b. 

In conclusion of this discussion, note that the formal- 
ism of developed pairing we use with particle number con- 
servation on average, as is known 22], works worse in the 
vicinity of magic numbers. For completeness, we include 
into analysis nuclei ^"^^Pb and ^^^Pb "dangerous" from 
this point of view. However, in this case it is natural to 
expect significant errors induced by this approximation. 

Let us go to odd-proton nuclei. Quadrupole moments 
of odd-proton neighbors of even tin isotopes, those of 
In {Z = 49) and Sb (Z=51) are displayed in Fig. E 
The Bogolyubov factor in ((TO)) is now reduced to -1 for 
the In chain and +1 for the Sb chain, thus, all differ- 
ence between predictions of the two functionals we use 
may come only from different values of the matrix ele- 
ments V\ of the effective field. We see that in the mid- 
dle part of both chains \Q\\ for the DF3 functional ex- 
ceeds that for the DF3-a in two times approximately. 
Remind that the effective field is that in the even Sn 
core, the same as was used for odd tin isotopes. Now 
the matrix elements are taken from the proton compo- 
nent Vp{r) whereas previously we dealt with the neutron 
component Vn{r). Both were displayed in Fig. [Hand 
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Te Xe Ba Ce Nd Sm Gd Dy 
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A 

FIG. 7: (Color online) Quadrupole moments of odd-neutron 
neighbors of even TV = 82 isotones with DF3 and DF3-a func- 
tionals. Triangles with the bars indicate experimental data 
i- 



TABLE II: Quadrupole moments Q (e b) of odd-neutron nu- 
clei in the state A. Theoretical values Qth and differences 
SQ = Qth — Qcxp are given for the functional DF3-a. 



nucleus 


A 




Qth 


SQ 


^"Ca 


1^3/2 


0.036(7) 


+0.040 


0.004 


«Ca 


0.040(6) 




0.000 


1/7/2 


-0.090(2) 
-0.066(2) 


-0.078 


0.012 
-0.012 


s^Kr 




-0.080(8) 




0.002 


159/2 


-1-0.443(3) 


+0.507 


0.064 


»^Kr 


2rf5/2 


-0.30(3) 


-0.355 


-0.06 


«^Sr 


lff9/2 


-H0.33(2) 


+0.335 


0.01 


«^Sr 


2d5/2 


-0.271(9) 


-0.245 


-0.026 


«9Zr 


I59/2 


+0.28(10) 


+0.262 


-0.02 


"Zr 


2rf5/2 


-0.176(3) 
(-)0.257(13) 


-0.195 


-0.019 
0.062 


i««Sn 




-0.206(10) 




0.011 


2rf5/2 


+0.31(10) 


+0.250 


-0.06 


"^Sn 


lff7/2 


+0.18(9) 


+0.029 


-0.13 


""Sn 


Iff 7*/ 2 


0.26(3) 


+0.377 


0.12 


"«Sn 


2^3/2 


0.094(11) 


-0.035 


-0.129 




-0.065(5) 
-0.061(3) 




0.030 
0.026 


i^^Sn 


2rf3/2 


-0.02(2) 


+0.063 


0.08 


""Xe 


2rf3/2 


+0.214(7) 


+0.217 


0.003 


»^Xe 


2/7/2 


-0.48(2) 


-0.376 


0.10 


"^Ba 


2d3/2 


+0.245(4) 


+0.254 


0.009 


"^Ba 


2/7/2 


-0.573(13) 


-0.445 


0.128 


i«Nd 


2rf3/2 


+0.32(13) 


+0.289 


-0.03 


i«Nd 


2/7/2 


-0.61(2) 


-0.518 


0.09 


"^Sm 


2rf3/2 


+0.4(2) 


+ 0.296 


0.1 


"^Sm 


2/7/2 


-0.60(7) 


-0.537 


0.06 




3p3/2 


-0.08(17) 


+0.195 


0.27 


Pb 


3p3/2 


+0.08(9) 


+0.272 


0.19 


201 T-ii 

^ Pb 


2/5/2 


-0.01(4) 


+0.137 


0.15 


203 DU 

Pb 


2/5/2 


+0.10(5) 


+0.284 


0.18 


205 pr, 
r^D 


^/5/2 


1 n OQ/'/i\ 


+U.00D 


n no 
u.uy 


209 pv 


2ff9/2 


n 

-U.o[Z ) 


-U.ZD4 


n 1 
U.l 


211pu 


2ff9/2 


+u.uy(D ) 


-U.zoo 


-U.o / 


1 I'j ^ 

"•"Sn 


i'*ri/2 


0.41(4) 


-0.776 


-0.37 




0.48(5) 




-0.30 


bn 




0.38(6) 


-0.703 


-0.32 


117o 

Sn 


1^*11/2 


-0.42(5) 


-0.593 


-0.17 


119c 

Sn 


1^*11/2 


0.21(2) 


-0.469 


-0.25 


121(3 

bn 




-0.14(3) 


-0.293 


-0.15 


123c 

Sn 


I/I1I/2 


+0.03(4) 


-0.123 


-0.15 




1^11/2 


4-0 1 (9) 


+0.039 


-0.1 


""Xe 


l'ill/2 


+0.62(2) 


+0.504 


0.12 


"^Ba 


l'ill/2 


+0.78(9) 


+0.588 


-0.19 


"^Gd 


1^13/2 
1^13/2 


-0.73(7) 


-0.791 


-0.06 


isipb 


+0.085(5) 


+0.0004 


-0.085 


193p,3 


-'-*13/2 


+0.195(10) 


+0.335 


0.140 




-'-*13/2 
1^13/2 


+0.306(15) 


+0.689 


0.383 


197p^ 


+0.38(2) 


+0.980 


0.60 


20"Pb 


-'-*13/2 


0.30(5) 


+0.665 


0.37 



we have seen that both of them are approximately two 
times bigger for the DF3 functional. This was caused 



8 



by a too strong collective 2^ state for this functional. 
For tin isotopes, this deficiency was partially compen- 
sated with different values of the Bogolyubov factor, see 
Tab. m Now we see the effect of the effective field di- 
rectly, without any distortions. In the case of the DF3-a 
functional, the rms deviation is worse than for tin iso- 



topes but moderately, {SQ)^^^[I)F3 - a] = 0.249 e b. 
As to the DF3 functional, the error grows drastically, 

\/('5Q)?ms[DF3] = 0.688 e b. There is a feeling that for 
the DF3-a functional the effective field Vp{r) is also too 
strong and the agreement will become better if it will be 
reduced by 20-30%. Within the mean field approach we 
use it can be achieved only with such variation of the 
central part of the EDF, common to DF3 and DF3-a, 
which changes the LM amplitude Tnp- 



a 



In 



— ■ 

.^/B » + , s ♦ 




■ 


9/2* DF3-a 


— □— 


9/2* DF3 


• 


7/2* DF3-a 




7/2* DF3 


• 


5/2* DF3-a 




5/2* DF3 




2.0 
1.5 
1.0 
* 0.5 
0.0 
-0.5 
-1.0 
-1.5 

101 105 109 113 117 121 125 129 133 

A 

FIG. 8: (Color online) Quadrupole moments of odd-proton 
neighbors of even tin isotopes with DF3 and DF3-a function- 
als. Triangles with the bars indicate experimental data J,]. 

Quadrupole moments of the odd neighbors of even lead 
isotopes are presented in Fig. [HI We see that here the 
difference between DF3 and DF3-a functionals is not big, 
especially in the right half of the figure where some ex- 
perimental data exist. Unfortunately, their number is 
limited, only one for Tl and five for Bi. In the last case, 
there are two sets of the data Both the calculations 
agree with one of them and disagree with the other. In 
such a situation, it is not worth to estimate the average 
disagreement numerically as it was made above. 

Let us go to the chain of odd-proton with magic neu- 
tron number N = 50, eleven nuclei from ^^Cu till ^^In, 
see Fig. [TUl Note that all of them, except ®^Y, are /3- 
unstable, many being far beyond the /3-stability valley. 
Here we deal with proton pairing, and the Bogolyubov 
factor in Eq. (jlOp varies strongly from one nucleus to 
another and depends on the state A. However, for this 
chain we obtain rather close predictions of both func- 
tionals used. In this case, the difference between both 
hmctionals is not as strong as in the case of the l/iii/2 
intruder state in the tin isotopes. Unfortunately, the ex- 




FIG. 9: (Color online) Quadrupole moments of odd-proton 
neighbors of even lead isotopes with DF3 and DF3-a func- 
tionals. Triangles with the bars indicate experimental data 




3/2 DF3-a 
3/2' DF3 
-•- 5/2' DF3-a 
-o- 5/2 DF3 
-♦- 9/2* DF3-a 
9/2* DF3 
Nb Tc Rh Ag In 



79 81 83 85 87 89 91 93 95 97 99 
A 



FIG. 10: (Color online) Quadrupole moments of odd-proton 
= 50 isotones with DF3 and DF3-a functionals. Triangles 
with the bars indicate experimental data [ij]. 



perimental quadrupole moment is known for one nucleus 
only, *''Rb. For this, both theoretical predictions practi- 
cally coincide with the experimental values. 

Quadrupole moments of the chain of odd-proton with 
magic neutron number = 82 are displayed in Fig. 
[TT] Here we deal with seven nuclei, from ^■^^I till ^^^Tb. 
Among them there are two stable ones, -'^^^La and ^^^Pr, 
the others are /3-unstable, but only ^'^^I lives hours, the 
rest, days or years. For this chain, four quadrupole mo- 
ments are known. For the DF3-a functional, disagree- 
ment with the data doesn't exceed 0.1 e b. For the 
DF3 one, disagreement is greater but not significantly. 
Comparison between predictions of the two functionals 
shows very close results for the 2^3/2 and 2^5/2 states. It 
demonstrates close results for the effective field Vp{r). A 
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FIG. 11: (Color online) Quadrupole moments of odd-proton 
N = 82 isotones with DF3 and DF3-a functionals. Triangles 
with the bars indicate experimental data m. 



stronger difference for the Ig^ /2 state is due to the differ- 
ence of the Bogolyubov factors in Eq. (|10l) similar to that 
we have seen for the neutron 1 177 /2 state, see Fig. [H It is 
worth to mention that the difference between two predic- 
tions for the proton intruder 1^13/2 state is significantly 
less than for the neutron one in the lead isotopes, Fig. 
[SJ More exactly, for neutrons the difference was strong 
for lead isotopes lighter of ^°°Pb and rather small, for 
heavier ones. Such behavior is explained with dynamics 
of variation of the spin-orbit density with filling the neu- 
tron shell which strongly influences the position of high 
j-levels. 

All quadrupole moments of odd-proton nuclei with 
known experimental values are collected in Table IIIII 
Just as for neutrons, we concentrated on the DF3-a func- 
tional which turned out to be more successful than the 
DF3 one. For the DF3 functional, the average error for 



protons is too big, J {SQ)'^^-^g[DFS] = 0.589 e b, whereas 



((5Q)2^,[DF3 - a] = 0.254 e b. However, for the DF3- 
a functional in average protons are also described worse 
than neutrons. The main contribution to this deviation 
comes from In and Sb isotopes, odd neighbors of even 
tin nuclei. As it was written above, it is the result of too 
strong quadrupole field Vn.p{r) even for DF3-a functional, 
not only for the DF3 one. For neutrons this drawback 
is partially hidden with multiplying by the Bogolyubov 
factor, see Table U but for protons it appears to the full 
extent. In many cases, we feel, the main reason of dis- 
agreement is the neglect of corrections from the particle- 
phonon coupling. Although in semi-magic nuclei this 
coupling is relatively small and can be accounted for with 
perturbation theory, sometimes these corrections could 
be noticeable. E.g. for the excited 3/2+ state in ^osrpj 
a strong admixture of the (2+ l/2"'')3/2 state may be 
expected. It could explain the very strong deviation SQ 
in this case. 



IV. CONCLUSIONS 

Quadrupole moments of odd neighbors of semi-magic 
lead and tin isotopes and iV = 50,iV = 82 isotones are 
calculated within the self-consistent TFFS on the base 
of the Energy Density Functional by Fayans et al. Two 
sets of the EDF parameters are used fixed previously. 
Namely, the functional DF3 [^,0| is used and its modifica- 
tion DF3-a nM ■ The latter is obtained with a small vari- 



TABLE III: Quadrupole moments Q (b) of odd-proton nuclei 
in the state A. Theoretical values Qth and differences 5Q = 
Qth " Qcxp are given for the functional DF3-a. 



nucl. 



X 



5Q(DF3-a) 



Id. 



'Sc If; 



3/2 



7/2 



'Rb 2p3/2 



lo^ln 
i°^In 
i^^In 
"^In 
"^In 



"^In 
"^In 
i^ifn 

i25ln 
i^^In 
"^Sb 
"^Sb 
i^^Sb 

123 



lS9/2 
lS9/2 
I59/2 
I39/2 
I59/2 
I59/2 

I59/2 
I59/2 
139/2 
I59/2 
I59/2 
lfl9/2 



5/2 
5/2 
5/2 



Sb lgr/2 
"'Cs lgr/2 
""La 137/2 



^Eu 2d. 
^Tl 3d: 



'Bi Ih 



'Bi Ih. 



5/2 
3/2 
9/2 

9/2 

■9/2 



^Bi lha 



9/2 



^Bi Ih 



9/2 



0.0585(6) 
-0.156(3) 
0.120(6) 
0.168(8) 
+0.134(1) 
-f0.138(l) 
-K0.83(5) 
+0.81(5) 
+0.84(3) 
+0.80(2) 
+0.80(4) 
+0.81(5) 
0.58(9) 
+0.829(10) 
+0.854(7) 
+0.814(11) 
+0.757(9) 
+0.71(4) 
+0.59(3) 
-0.36(6) 
-0.37(6) 
-0.36(4) 
-0.45(3) 
-0.49(5) 
+0.051(1) 
+0.20(1) 
-0.077(6) 
-0.059(4) 
+0.29(2) 
+0.74(15) 
-0.93(7) 
-0.68(6) 
-0.81(3) 
-0.59(4) 
-0.76(2) 
-0.55(4) 
-0.60(11) 
-0.516(15) 
-0.37(3) 
-0.55(1) 
-0.77(1) 
-0.40(5) 
-0.39(3) 
-0.83(5) 
-0.60(5) 



0.069 
-0.139 



+0.132 

+0.833 
+0.976 
+1.113 
+1.165 
+1.117 
+1.034 

+0.965 
+0.909 
+0.833 
+0.743 
+0.663 
+0.550 
-0.882 
-0.766 
-0.721 

-0.739 
-0.031 
+0.103 
-0.120 

+0.156 
+0.227 
-1.323 

-0.945 

-0.454 



-0.342 



-0.508 



0.010 
0.017 
-0.019 
0.029 
-0.002 
-0.006 
0.00 
0.17 
0.27 
0.36 
0.32 
0.22 
0.45 
0.136 
0.055 
0.019 
-0.014 
-0.05 
-0.04 
-0.52 
-0.40 
-0.36 
-0.27 
-0.25 
-0.080 
-0.10 
-0.043 
-0.061 
-0.13 
-0.51 
-0.39 
-0.64 
-0.14 
-0.36 
0.31 
0.10 
0.15 
0.18 
0.03 
0.21 
0.43 
0.06 
0.05 
0.32 
0.09 



10 



ation in the initial functional DF3 of the spin-orbit terms 
and rather strong, of the effective tensor term. The latter 
influences the spin-orbit splitting in nuclei with partially 
filled spin-orbit doublets where the spin-orbit density psi 
is non-zero. Due to different spin-orbit splitting, the po- 
sition of high j-levels and all the single-particle spectrum 
in the vicinity of Fermi level are often different for two 
functionals under consideration. The complete set of the 
QRPA-like TFFS equations for the effective field is solved 
in a self-consistent way. 

Recently, we used this method (T^ to calculate charac- 
teristics of the 2^ states in even tin and lead isotopes with 
the use of the DF3-a functional. In this work, quadrupole 
moments were calculated as well of odd tin and lead iso- 
topes and odd-proton neighbors. In and Sb isotopes in 
the first case and Tl and Bi isotopes in the second one. 
Strong sensitivity of quadrupole moments of odd- neutron 
tin and lead isotopes to the single-particle spectrum in 
the vicinity of the Fermi surface was found in 10] due to 
the Bogolyubov factor u\~ vf — {e\ — ^)/E\. 

To examine this effect in more detail, in the present 
paper we calculated quadrupole moments with the DF3 
functional. These calculations confirmed the observa- 
tion of Ref. |10[. In addition, it was found that the 
static quadrupole effective field V{r), the main ingredi- 
ent of the formula for the quadrupole moment Q\, also 
can be significantly different for these two functionals, 
again due to different level structure in the vicinity of 
the Fermi level. This occurs in the middle part of the 
tin isotopes where the DF3 effective field approximately 
two times stronger than the DF3-a one. The latter looks 
more realistic as follows from the analysis of the charac- 
teristics u}2 and BE2 in even tin isotopes. [23!] For odd tin 
isotopes, a too strong effective field is hidden partially 
with smaller values of the Bogolyubov factor, but the 
situation is worse for odd-proton neighbors. In and Sb 
isotopes, where the DF3 functional leads to a strong dis- 
agreement with experiment. For other chains considered 



the difference between the predictions of both functionals 
is less, but on the whole, the DF3-a functional describes 
quadrupole moments better than the DF3 one. 

For the DF3-a functional, the agreement, on average, 
can be considered as reasonable. For 42 quadrupole 
moments of odd-neutron nuclei, the average disagree- 
ment between theory and experiment is not so small, 

\J (5Q)rnis — 0.189 e b. However, it is concen- 
trated mainly in 15 intruder states for which we have 

\/(^<3)?ms [intruder] = 0.269 e b. For the rest of 27 
"normal" states, the disagreement is rather moderate 

\/(^<3)?ms [normal] = 0.125 e b. 

For protons, agreement is worse. Leaving aside spe- 
cific reasons in several cases, we see the main source of 
disagreement here, just as for neutrons, in neglecting the 
effects of the coupling of siiigle-particle degrees of free- 
dom with phonons, see e.g. [1, UHII- 

The comparison between two functionals favors DF3- 
a which has a rather strong effective tensor force. For 
further improvements, phonon contributions should be 
taken into account which would pave the way to extend 
the present investigation to nuclei with both non-magic 
proton and neutron subsystems. 
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